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FOREWORD 


This report, which presents working rules and recommendations for flutter 
analysis using the aerodynamic-influence-coefticiant method, was prepared by 
the Aeroelastic and Structures Research Laboratory, Massachusetts Institute of 
Technology, Cambridge 3¥, Massachusetts tor the Aircratt Laboratory, Wright 
Air Levelopment Center, Wright-Patterson Air Force Base, Ohio. The work was 
pertormed at the MIT unde: the direction of Professor H Ashley, and the 
Project wag supervised by Mr. G, Zartarian. The research and development work 
was acconplished under Air Force Contract, No, AF 3S (816) ~2uae, Project No, 

1970 (Unclassi tied Title) “Meroelastici ty, Vibration and Noise," and Task No. 
134 Nclassified Tittle) "Theoretical supersonic Flutter Studies, * Mr, 

Walter J. Mykytow of the Dynamics Branch, Alroratt Laboratory, is task engineer. 
Research was started on | July i984, This report is part of a continulng effort 
in the flutter of aircratt structures at supersonic speeds. This is Part {I 

of this report which I$ published In two Separate parts, 


The author, Mr. G, Zartarlan, Is indebted to Protessor H, Ashley, Or, P, T, 
Hsu and Mr, A, Heller for thelr contributions to the research. In addition, 
‘acknowledgements are due to Mr. G. Anitole for Preparing the figures and to 
Mrs. B. Marks for typing the finai manuscript, 


This document, including the Illustrations, is class! fied CONFIDENTIAL 
(excepting the title) because It contains the development of improved methods 
tor conducting Supersonic flutter analysis; hence more accurate flutter 
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ABSTRACT 


On the basis of investigations described in » preceding re- 
port (Ref. 1) and subsequent research, recommendations are thade 
concerning methods of application of aerodynamic pressure or 
velocity potential influence coefficients to practical flutter 
analyses. Following a brief discussion of the relative merits 
of three basic types of elementary areas, working rules are 
given for ae Plenforms of principal interest. Simple illua- 
trative examples have been included to clarify some of the de- 


areas, which is deemed to be the most sat cain gay from the over- 
all etandpoint, formulas and tables of certain relé 
are presented which will facilitate future tabulations of aero- 
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SECTION I 
INTRODUCTION 


In a preceding report (Ref. 1), the need for direct, sye- 
tematic numerical methods for the determination of supersonic 
tor aaae aerodynamic forces on wings was discussed with partics- 
lar emphasis on the case of simple harmonic motion, which oc- 
curs in flutter calculations. Such an epproach involves a divi- 
sion of the planform into Small areas and use ef so-called sero 
dynamic influence coefficients. Threa alternative types of 
elementary area were investigated with the ultimate aim of find- 
ing the Particular one which would yield the most roe gd a 

rom an overall} Standpoint. Some of the @spects of this 
problem which had to be considered were the adaptab lity of the 
various grids to Planforms of interest; limitat ong on the ele- 
mentary box size, reduced frequency of tke oscillations and Mach 
tumber; tha relative @ase of Cabulating the aerodynamic coef - 
ficients; and finally the accuracies attainable in the general - 
ized forces which enter the flutter equations. 


In the assumed mode “pproach to the flutter Problem, the 
uantities directly involved are the aforementioned 8eneralized 
orces and not the pressure distributions themselves. There- 

fore, the ultimate aim is to determine the generalized forces 
satisfactorily. Evidently, if the pressure distributions are 
found accurate] everywhere on the Planform, the resulting gene- 
ralized forces which are weighted integrals over the planform 
area of these pressure distri utions) will also be accurate. 
Cases ape cwever, when the pressures at various control 


rpereciably from the true pressures. Fortunately, because of 
the averaging effects of numerical inte rations over the plan- 
form, the resultant generalized forces in such Cassese will still 
be sufficiently accurate, provided the recommended rules in Sec- 
tions III and Iv are followed, 


In Ref. 1 the primary interest was centered on the appil- 
cability of the so-called aerodynamic pressure influence co- 
efficients (abbreviated from here on as PIC), becaugce this type 
of influence coeff‘ etent appeared to be the most straight forward 
in applications. After completion of the ressarch reported in 
Ref. 1, it was found that to determ’ ne generalized forces a 
better approach would be to make use of the velocity potential 
influence coeffictent (abbreviated from here on as VIC). Im ell 
the applications made up to date, the VIC-method has prover it- 
self to be superior to the PIC-method, and has indicated the 


Larne nT 
Manuscript released by the author February 1956 for publication as a WADC 
Technical Report. 
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following two important adventages: 


(1) The Vic-tabulation Would require less than half the 
machine operat ions than those for the PIC. 


(2) For planforme with subsonte edges, and for a given 
accuracy, fewer boxes ere heeded with the ViC~method, 


Because of the importance of thig new @pproach, it is reported 
here in advance, pending the vere atten of detailed investi - 
ations on the ViC-method. Further applicationa of thig method 
will be undertaken and reported on in the future. Although the 
rules stated in Sections Iil and tv are based riwarily on the 
findings about the PIC-approach, they are still 

modifications when employed in connection with the VIC 


Inasmuch as the so-called Mach box grid system offers ad- 
vantages in many reg ects, especially for Planforms with gub- 
sonic edges at M< V5 » the emphasis has been Placed on the use 
of this system. Kules have been devised to assure setisfactory 
results, with the ma X Limam gpg ed and Simplicity, whenever 
possible. Additional rules and techn 
to meet the requirements cf specific plantorms and edge con- 
figurations. ection IV 
to account for the effect of the singular downwash at side edges 
on the airloads. Such refinements, which improve markedly the 
accuracies of the final results (Ref, 1), have been recommended 
Vor use Whenever they are needed and can be incorporated con- 
veniently, 


Although not specifically stated, some of these rules are 
applicable to the Square grid system, for which tables of the 
Pic's are already available (Ref. 2). 


At relatively high Mach numbers, say 52.5, the recently 
developed apiscon theory" is Promising for many types of plan- 
forms. Preliminar investigations (Refs. 3,4) have indicated 
that, at these Mac numbers, airloads on two-dimensional or 
purely supersonic airfoils with small to moderate sweep can be 
predicted "ha. An extension to the three-dimensional 
case seems featible. The accuracy of such predictions improves 
as the Mach number ig increased. In view of the simplicity of 
this theory, its use, whenever applicable, ia preferred for this 
high Mach number range. 


In Appendix D, a simple flutter problem is formlated to 
illustrate the various steps in applying either of the two types 
of influence coefficients. If the £ utter probiem is to be 
solved by means of a high-speed digital computer, two alterna- 
tive ways can be followed in feeding the aerodynamic influence 
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coefficients inte the wachine: (1) to store the required table 
of influence coefficients which have already been calculated 
elsewhere, or (2) to store the sub-routines which generate the 
influence coaf icients whenever they are called for. The final 
cholee as te which precesa should be coniages depends mainly on 
the storage Capacity and the speed of the Computer, and the aage 
of PCOgTamming . 


This report is intended primarily to supply veck ing rules 
in the use of the Mach rid system for aerce asticieus interes- 
ted in solving practice flutter problems. Ko effort has been 
made here to gu stantlate these rales, since they have been ex- 
tensively discussed in Ref. 1. It should be emphasized here 
that these rules were arrived at by experience on @ limited num- 
ber of representative trial cases, with due Consideretion being 


of the method, Therefore, these rulee do not guarantee a given 
accuracy, and their use does not necessarily result in « unifors 
accuracy for all types of Planforms. "hey are expected, however, 
to yield acceptable calculated airload for flutter analysea. 


With the VIC-method which is concluded to be the better ap- 
proach, it is felt that ft is now unnecessary to tabulate aero~ 
dynamic influence coefficients, rkcn one of the tabies in this 
report (page 45), the procedure for obtaining the VIC's according 
to the approximate expressions is quite stra ght forward, Engi- 
mers using the present method can have these tabulated in a very 
short time for any Mach number and reduced frequency. Most 
flutter calculations will be carried sut on high-speed digital 
Computing machines, on which very simple subprograms can be de- 
vised for instantaneous Senerations of the coef iclents as needed, 


Under an extension of the subject contract, the following in- 
vestigations are contemplated: 


(1) The possibility of using a subdivision technique for partial 
boxes at the leading and trailing ed-'s, 


(2) A simplified procedure for handliny the downwash singularity 


ment of this type for subsonic leading and side edges. 


{3) Applications of these methods to wings of various planforms 
Straight tapered, delta and swept), for which experimental data 
are available, ‘ 
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SECTION If 
SUMMARY OF THE RELATIVE MERITS OF THE 
THREE BASIC GRiD SYSTEMS 


In a preceding report (Ref. 1), several numerical techniques 
for the ca culation of airloads on oscillating wings in super - 

sonic flow Were atudied, all of which have their et in a 

suggestion by Pines (Ref, 5), Briefly, this approac starts 

with the definition of the aerod mic influence coefficient, 

according to which the Pressure (or the velocity Potential) at 

any point (x,% ) on tne planform is related to the values of the 

vertical velocity (“downwash"') on elementary areas ("boxes") lo- 

catec in the region bounded by che forward Mach lines emanating 

from (x, % ). With certain apecific exceptions, auch as the down- ‘ 
wash Singularity iuert off a side edge, each aerodynamic influence ‘ 
coefficient represents the pressure (or the velocity potential) 

induced at (~~) by a Constant unit downwash over a Particular 

elementary area, with zero downwash over all other guch areas. 

In view of the linearity of the theory, the Cotal pressure (or 

the total velocity Potential) at x4 5 

to the sum of the influence coeff cients for all boxes affecting 
that point, each one weighted with an ’ppropriately averaged 
downwash for its box. en the quantity of interest at (x,% ) 
ie the pressure, one uses the preasure influence coefficients 
(PIC). ‘When this quantity is the velocity potential (¥Y ) one 
must employ the velocity potential influence coefficient (vic). 


Three basic shapes of elementary areas were investigated, 
in connection with the PIC-method, with the objective of deter- 
mining the overal] practicabilicy of each in numerical applica- 
tions, 

1. The square box (the shape originally proposed by 

ines); 
@. The Mach box, which is @ rectangle with its 
diagonals paral)! Co the Mach lines; 


3. The characteriat Le box, which is a rhombus 
bounded by four Mach lines. 


These alternative schemes were applied to simple trial cases 
for which analytic solutions are known. In go doing, several 
relevant points were established regarding the Suitability of 
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each of the grid systems. seyentagne and disadvantages @8£0c~ 
lated with each one were recognized. 


Before et Seupeing to apply the influence coefficient methods 
to particular Problems, it is desirable to underetand Clear] 

the relative Werite of each box shape, so as to make a Judielous 
choice, For this purpose a Summary of the conclusions reached 
according to the investigations for the Pressure influence coef. 
ficiente (Ref. | is now Presented, In, Beneral, the conclusions 
reached there ere expected to Spply for the VIC-method as well, 


II.l The Square Box 
Advantages: 


tics Tabulations 5« the preaaure influence coefficients 
are alresay available for a 26 x 20 grid eyuten (Ref. 2), 
The number of boxes that can usually be taken in the 
chordwise direction (20) is large enough for most pur- 

- These boxes may nearly always be used satisfac - 
torlly for Purely supersonic planforms (e-g., the wide 

a They may also be used for Planforms with aub- 
Sonle edges when M>J/2. 


2. The grid System, and therefore the pressure (or 
velocity Potential) and downwash Control points, re- 
main fixed when the Mach number is varied. This is a 
desirable feature Since it facilitates the necessary 
numerical integrations for the Beneralized forces. 


Disadvantages: 


l. This System cannot be used for Planforms having sub- 
Sonic edges when Ma«vV/z. 


2. On the basis of trial “xamples with the VIC, this 
method ig very likely to SUupersede entirely the origi- 


2: The effect of the downwash singularity near a side 
edge cannot be easily incorporated, when this refine-~ 
ment is needed. 


s, Tf exact expressions (as far as the frequency is 
concerned) for the influence coefficients are uged in- 


Stead of the approximate ones giver, by Pines (Ref, 6), 
fewer boxes may be sufficient n some instances. 
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Tabulations using these exact expressions are much more 
complicated, however, than the corresponding ones for 
the Mach box (see Appendix A), 


5. Unlese one tabulates the contributions due to "half 
square boxes," the condition of sera pressure (or equi- 


On @ Streamise side edge, with the result that lift 


6. The centers of boxes do not fall on the Mach line 
which separates the purely supersonic and the mixed 
regions near a side edge. A refined chordwise inte- 
gration technique for calculation of Beneralized for- 
ces (with the PIC-wethod), as illustrated in Appendix 
C, therefore cannot be used, 


Il.2 The Mach Box 
Advantage : 


i. This system can be used conveniently and with satis- 
factory accuracy for planforme having subsonic edges 
and for all supersonic Mach numbers, (4>a>s2 )., 


2. The singularity in downvash at e side edge may be 
accounted for. 


3. Tables based on Watkin's exact expressions for the 
influence coefficients (PIC or VIC) are rather simple 
to compute (See Appendix A). 


4. Pressures are obtained on the Mach line which sepa- 
rates the purely supersonic and the mixed regions near 
a eide edge. This is desirable if the inte ration 

technique of Appendix ¢ is needed (for the IC-method) . 


Disadvantages: 


1. Unless one tabulates the contributions due to "half 
Mach boxes," the condition of sero pressure (or zero 
velocity potential) cannot be satisfied on a side edge, 
with the result that lift and moment distributions can- 
not be made to vanish identically at auch an edge. 


2. The control point locations vary when the Mach num- 
ber is tre ed and hence the downwashes must be re- 
~alculated for each Mach number. However, this task is 
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vised for ®ach Mach num er; again, thie ig relatively 


Simple and Constitutes only a very small Percentage of 
the total computational efforc, 


boxes need NOt be too large, Also, for MP25 , it 
pecoumended that the much simpler "piston theory" 
8, 


3,4) be used whenever “pplicable, 


Il.3 The Characteristic Box 
Advantages; 


lL. Thig System can be used for Planforms having Sub- 
Sonic ee for ail Supersonic Mach numbers, 
Fins.2 ), 


2. The singularity in downwash SC @ aide edge may be 
@Ccounted for. 


OF Tabulations (PIC or VIC) are not *xpected to be tog 
diffierrt , although they will be more labor ioug than 
those for Mach boxes, 


4. The Condition of =*ro pressure (or zero velocity 
Potential) can be Satisfied at discrete Points on a 
8ide adge. 

5. Pressures @re obtained on the Mach line Separatin 


the Purely eu ersonic and the mixed regions near a bide 
edge (PIC-met od). 


Disadvantages; 


1. The Control poinr locations “yr A when the Mach num- 
ber ig changed. (This is » minor o jection ag Stated 
in connection with the Mach box, ) 


2. As in the case of the Mach box, for Mogg 


&panwise boxes must be taken to ensure a #ufficienr num- 
er of boxes in the Chordwige direction, 


CONFIDENTIAL 


CONFIDENTIAL 


3. Although ac first glance this box system might ap- 
pear quite versatile for awept leading edges, experience 
has shows that ( exeepe for & leading edge exactly paral- 
lel to a@ Mach line) it yields poorer results than when 
the leading edge is represented by a jagged boundary 
composed of Lines Bespseeten tor and parallel to the flow. 
To represent the lead ng edge in the latter fashion, 

when characteriatic boxes are i eo employed, tables 
must be calculated of aerodynamic influence coefficients 
for "half-rhowbic" or triangular boxes. Furthermore, 
experience indicates that when a bet has a strearwise 
side edge or a swept leading edge and the downwash varies 
appreciably in the chordwise direction, it is desirable 
to have influence coefficients tabulated for "quarter- 
rhombic’ boxes. These requirements increase the total 
number of entries in the tables by a factor of four, 
which is a serious drawback relative to the other two- 
pee sine These problems are discussed in detail in 
Ref. 1. 


It should be noted that some of the above advantages and 
disadvantages are based on che condition of keeping to a minimum 
the computational labor to determine satisfactorily the gene- 
relized forces; (this implies keeping the number of boxes on the 
planform to a minimum), If many more than the required number 
of boxes are taken, some of the refinements, such as the inclu- 
sion of the effect of aide-edge downwash Singularity or the use 
of the more elaborate Stages ies Cechniques of Appendix C, may 
become unnecessary. There ore, for theze cases, some of the ad- 
vantages and disadvantages mentioned above may become inaignif i - 
cant. 

From experience to date, it Sppears that the Mach box system 
rates most favorably for the majority of practical applications, 


Accordingly it is recommended as the basia for future tabulations 
cf supersonic aerodynamic influence cesffictente, 
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SECTION IIt 
RULES FOR APPLICATION OF THE MACH BOX SYSTEM TO 
PLANFORMS OF PRACTICAL INTEREST 


In order to carry out flutter aps hon efficiently using 
the Mach box scheme, ¢ ‘rtain working cules. which are given be- 


low for specific Planforms and edge conditions, must be followed. 


It is possible to enunciate more genere) rules which apply for 
all combinations of planform and Mach number. The latte: ret 
would allow the analyst less freedom of choice in a par’ leular 
case, Sowaver, since it would overlc ok some of the possible eim- 
plifications and refinements, 


Planforws of current interest may be classified into three 
categories: the straight tapered, the sweptback tapered, and 
the triangular. As che Mach number ig varied, each of these 
plentorne may exhibit supersonic or subsonic edges, the former 

eing characterized by a sweep angle less than that of a Mach 
line, and conversely. For different ranges of Mach number, it 
ig necessary to follow different sets of rules. 


Although these rulee are based on the ex erience gained by 
carrying out examples using the PIC's (Ref. iss a few trial cai- 
culations on some of these examples with the use of the VIC's 
have indicated that most of the rules still apply, with increas- 
ing resultant accuraciles for planforme with subsonic edges. 
These rulee should therefore be adopted for both methods, unlegs 
specifically stated otherwise. 


IIl.1 Purely Supersonic Pianftorms 


The following rules apply for purely supersonic planforme 
(planforms with no subsonic eading or trailing edges and no 
side edges; e.g., the wide delta). 


Rules: 


1. The dimensions of the elementary areas to be distribu- 
ted over the planform are found from the condition that, 
for the highest Mech number to be considered in the analy- 
Siz, a minimum of eight boxes chordwise along the root 
chord is recommended. When peastene: the Bee dimen - 
sion of the boxes thould be kept constant, because then the 
spanwise locations of the control points do not vary, and 
the overall width of the grid of doxee bears a fixed 
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relation to wingspan. As the Mach number ig decreased, the 
resulting reduction of the chordwise box dimenaion & causes 
only a moderate variation of the reduced frequency 


his (wh jv)(o*/a*) 


for a given frequency of oaciilation ew . However, if the 
Fange between the maximum Mach number, ,., and the Mach 


number when any of the edges becomes sonic, Me, 1é large, 


the number of chordwike boxes may become quite high, If 
scent boxes are taken for Mmax , then for Meon there will 
e 


3 Vi -7 d's My? =/ boxes, 


If thie number becomes i tactical on account of Limita- 
Cions ef tabulatione or the capacity of the ¢ uting 
machine, @ reduction of the number of Spanwise boxes is 
recommended at the lower Mach numbers. in any Case, no 
leas than eight boxes along the root chord should be em- 
ployed. For most of these types of Planform the ht Sea 
tion of this rule will also ensure a sufficient number of 
Spanwise boxes, 


2, The number of chordwise boxes is further governed by 
two factors: the accuracy of the tabulations and the most 
sinuous assumed mode shape. 


a. If approximate formulas for the influence coeffic- 
dents (Appendix A) are used, the huuber of chordwise boxes 
must be chosen such that 


ie wh M* 
ae eee 0.3 
4 (ob 3 < Eq. (3.1) 


b. In the evaluation of the generalized forces, ag- 
sociated with the assumed mode type of flurter analysis, 
the mode which is the most sinuous in the chordwise direc- 
tion will determine the maximum tolerable box size. A mini- 
mum of six boxea per "“haif wave length" of this mode shape 
should be used. 


3. After the box size is determined ay conditione 1 and 2, 
@ grid system can be constructed. Altho 

for proper positioning of this etid system on the Planform 
cén be given for all cases, Che example shown in Fig. III.1 
will serve as a wide. The he edge of the actual plan- 
form will be replaced by a jagged line. Therefore, the 
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Fig, I. 1 Positioning of the Mach Box Grid System on g Purely 
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chordwise Positioning of the box system with respect to 

the leading edge must be auch thac the leading edge area 

of the planform left out by the grid System is a proximately 
equai to the aree outside the planform taken in y the grid. 
Since the trailing edge is Supersonic, no particular re a- 
tionship between chat line and the box centers adjacent to 
it need be imposed. 


Supersonic Planftorm 


4, According to Rule 3, the actual planform is replaced by 
an approximate one composed of boxes with Centers on the 
Planform. Near the Craiiing edge, all boxes partialiy on 
the planform, regardless of the locations of their centers, 
ere to be included, Two types of such partial boxes ap- 
pear. Referring to Fig. III.1, box @ hag its center on 
the planform, while box & hae ite center behind the trail- 
ing edge. For box a , the pressures (or the velocity 
potential 4 ) and the associated weighting factors for 
the generalized forces (Appendix C) can be computed Like 
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way 


those for any other complete box. for box & , the pres- 
sures (or 4y) and the weighting factors should be aseumed 
equal to the "Gael erent ones for the adjacent upstream 
box © . For each partial box, its contri ution to the 
seneralized force must be multiplied by the ratio of its 
srea on the planform to ite total area. In addition, when 
dealing with the ViC-method, it is necessary to determine 

AY at points on the trailing edge, This 1s accomplished 
by extrapolation. For instance, the Ay at point / ig 
found by linear extrapolation using the valuse st the cen- 
ters of boxes 2 and a@ ; the OY at point 2 ie found b 
Linear extrapolacion are the value of zero for point j 
and the calculated value for the center of box e , 


5. With the pattern of boxes over the pisofore completely 
established, the (complex) downwash amp 
of each box ig found from the mode shape of Simple harmonic 
motion. Appendix B discusses this operation. With the 
downwash amplitude at the center of each box known, the 
Complex) amplitude of pressure uifference or of the velo- 
elty potential difference between tt. Upper and lower eur- 
& 


(P= bb hh » P= 5-H) 
can be ie ee at each control point from either Eq.(A.6) 
or Eq. (A.8) of Appendix A, 


6, The rectangular integration rule is to be applied for 
Beneralized force integrations. The Feséures or ays and 
the wei hting factors are assumed to be constant over each 
box. The feneralized force can be computed from Eq. s6. 10) 
or Eq. (C.14) of Appendix C, or from equivalent formulas. 


Ill. Planforms with Supersonic Leading and Trailing Edges and 
with Side Edges Parallel to the Fright Diresrtoe ie 


Aasociated with any subsonic edge, such a6 a side edge, 
there are two Special regions to be considered, These are the 
#O-called "nixed region," which is the portion of the planform 
influenced by the subsonic edge and the 80-called “diaphragm 
region" off the planform where disturbances exiet. The condi- 
tion of sere Pressure diffecence must be imposed over the dia- 
phragm region, Near « side edge paralisl co the flow, this 
condition is identical with the requirement that the difference 


in the velocity potential be zero at al) points of the atepbragm. 
6 


From aither of ¢ 
region can be calculated, Furthermore, it is known that the 
downwash distribution has a square-root singularity as the sub- 
sonic edge in approached from the diaphragm region. Two 
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alternatives in calculating the effect of this diaphragm region 
on the mixed region may be followed: 


Alternative 1: If the unknown downwash over each eteenr die 
Ox 18 aesumed constant, the downvash over such boxes 

may be computed by satisfying the condition of zero 
pressure (or rero 4¢ ) at the centers. Once these 
downvashes are determined, the er Ha 0 boxes may then 
be treated like any other box on the Planform. As a 
Beneral rule, this is satisfactory when there is a mini- 
mum of ten boxes along the side edge. It should be 
pointed out that with ten boxes and the VIC-method, one 
obtains better accuracy than with the same number of 
boxes and the PIC-method, so that thie minimus may be 
lowered somewhat when the former is eaployed. 


Alternative 2: When fewer than thie minimum are employed along 
@ side, the assumption of conetant downwash over each 
diaphragm box adjacent to the side edge is not adequate 
unless the mixed region is small compared to the plan 
area and contributes a amall portion to the generalized 
force. In most practical cases, however, the ac gg 
near the side edge will be quite large since the maxi- 
mam deflections uring flutter occur near the tip. Then 
one must incorporate the effect of the downwash singu- 
larity. This procedure ig illustrated in Section IV. 


Rules: The rules given for rigs | supersonic planforme apply 
in general also to this class of planforms, These are 
supplemented by the following additional rules: 


L. The approximate box size dictated by the conditions of 
Section I1I.1 should be adjusted so that there will result 
an integral number of boxea across the span. 


@. Referring to Fig. 111.2, the actual side edge 2% in re- 
placed by a’#. Cons gugttly. the actual mixed region e& 
must be replaced by a‘#c’. 


3. The diaphragm region fad is to be represented by com- 
plete boxes if the downwashes over such boxes are assumed 
constant. Its one edge if no longer a Mach line but 
ses jagged line formed by the boundar 8 of che foremost 
OXeR, 


4, The pressures or the velocity potentials over the mixed 
region 2’#c’ are to be found according to Section Iv. 
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Fig. W.2 Positioning of the Mach Box 


with Supersonic Leading and Trailing Edges and with 
Side Edges 


>» If few chordwise boxes are 


DIAPHRAGM 
REGION d 


Grid System on g Planform 


employed near the side edge, 


the integration technique of Appendix C is to be used whe: 


Computing generalized orces, i 
tangular rule. Yhis refinement 
PIC-method, 


IT1.3 Planform with Subsonic Leadin 


When a planform with subsonic 1 
two serious difficulties arise. Fir 
leading edge by a jagged line causes 
pressure or of the velocity potentia 
each box, Secondly, in the di 


that near a side edge. Unlike the 


cannot be incorporated phar i a 
difficulties dictate the use of ala 


WARDS TR 56-97, Part 2 ~14- 


n place of the simple rec- 
1a necessary only with the 


and Supersonic Trailing Edges 
ae edges is considered, 
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aphragm i ahead of the leading 
edge, there exists a Singular downwagh d 


stribution similar to 
ide edge case, this effect 

in the thecry. Both of these 
Fge number of boxes, 
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Rules: 


1, A minimum of twelve boxes along the mid-span chord igs 
required. 


2. Rules 2, 3, & of Section ilI.1 apply. aepure IIL,3 ie 
furnished as an example of how a typical planform mighe be 
overlaid with Mach boxes, 


3. If side edges exist, the approximate box size dictated 

by previous conditions should be adjusted so that there 

Will result an integral number of boxes across the span. 

Again the spanwise imension should be Kept constant ae the 
ch numbers more than 

twelve chordwise boxes will then be Present, and this is 

desirable because the leading edge becomes "more subsonic." 


4, The diephragm region is comprised of complete boxes in 
ne eee region off the Planform, as shown in Fig. 


5. With the pattern of boxes over the plenfors completely 
established, the (complex) downwash amplitude at the center 
of each box on the planform is found from the mode shape of 
simple harmonic motion (Appendix B). 


6. With the downwash amplitude at the center of each box 
on the plenform known, the constant downwash amplitude over 
each box in the diaphragm region is determined from the 
condition of zero pressure or zero 4%) at the centers 

of diaphragm boxes, according tu Section Iv. 1. 


7. With the downwash amplitudes at the centers of all 
boxes known, the (complex) amplitude of pressure difference 
(or Ay ) between the Upper and iower surfaces can be com- 
puted et each control pours on the planform from Eq. (A.6) 
or Eq. (A.8) of Appendix A. 


8. The potentials at trailing edge points are found by 
extrapolation as in cule 4, Section IiI.1. 


9. The rectangular cule is to be applied for generalized 
corce integration, the pressures (or Ay ) and the weight- 
ing factora being assumed constant over each box. The 
enerelized forces can be Computed from Eq. (C.10) or 
q. (C.14) of Appendix c. 
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Fig.I.3 Positioning of the Mach Box 
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Grid System on a Plantorm 
with Subsonic Leading Edges and Side Edges 
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The foregoing rules were arrived ac originally in connec - 
tion with the PIC-method (Ref. 1). They were stated after giv- 
ing due consideration to the following factors: (1) accuracy 
for the generalized forces, (2) the amount of computational 
labor. Neediess to say, higher accuracy requires a larger amount 
of work. Therefore, for the final decision a compromise had to 
be made, i.e., one had to determine the point beyond which the 
improvement in accurecy ig not sufficient to warrant the assoc- 
lated inecreaees in the computations. 


With the same number of boxes, orsones has indicated 
that the ViC-method yields ei nificantly better accuracy then 
the PIC-method. In view of t is, one would think of reducing 
the minimum number of boxas set forth in Rule 1 from twelve to, 
say, eight. However, it is regardad as preferable to keep the 
miniwum number of tweive for the following reason. For this 
clase of planforme, the accuracies attainable for the general- 
ized forces are, in general, poorer than those for purely super - 
sonic planforma when the above rules are employed, By keeping 
the miniwum number of boxes at twelve with the ViC-method, the 
computational precision is brought more in line with that for 
purely supersonic planforms. 


TIl.4 Planforma with Supersonic or Subsonic Leading Edgee and 
Subsonic Tr ng Edges 


For planforms with subsonic frailing edges, the rules etated 
previously in Section III. 3 @pply in general with the following 
modifications to account for the diaphragm region behind the 
treiling edges (see Fig. [I1.4). 


Rules: 


1. If the leadin edge is supersonic, a minimum of twelve 
boxes along the m depan chord is to be used; if aubsonic, 
sixteen boxes. 


2. The diaphragm region behind the trailing edge must be 
composed of complete boxes only, Consequently the trailing 
edge is replaced uy &@ jagged line { see ig. ILI.4) and in 
th’s case partial boxes are not allowed. 


3. For diaphragm regions associated with a subsonic iead- 
ing edge or a side edge Parallel to the flow, the condition 
of 4¥=O0 is equivalent to the condition of sero pressure 
difference. For the diaphragm boxes behind a subsonic trail- 
ing edge, the equivalence of these two alternative condi- 
tions is no longer true. Rather than setting 4Y=0 for 
such diaphragm boxes, one wust proceed in the following 
manner to satisfy the true boundary condition of aph=oO , 
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Consider the center =f box 2 . The value of 4¥ at this 
point should be taken to be 
-ch (4a-%) 
(ay), = (ag), ee") Eq. (3.2) 


where (49) ie the potential difference at che upstream trail- 
ing edge point / and (%~% ) represents the distance be- 
tween the center of box @ and point / .# 


“It should be noted that (9), must be computed for the gene- 


ralized force calculations (Appendix C), so that this necessary 
step does not add much to the complexity of the problem. 
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SECTION IV 


DETERMINATION OF DOWNWASH DISTRIBUTION ON THE 


DIAPHRAGH REGION AND ITS EFFECT ON 
THE AIRLOAD DISTRIBUTION 


Aesociated with any bubsonic edge, there exists a disturbed 
region off the planform which influences the pressure distribu- 


tion on the planform, Following Evvard's concept (Ret. T), @ 
erweable diaphragm may be introduced in th 

This diaphragm is assume te coincide with a atream sheet, and 
therefore it will not alter the flow. It can sustain no pres- 
sure differance at any point between its Sop and bottom sur- 
The latter cordition allows the determination of the 
unknown vertical velocity (downwash) distribution on the dia- 


thin iwp 


facas, 


a new surface which is purel supersonic and for which the down- 
known everywhere. The pressures (or the velocity poten- 


wash is 


® region. 


The combination of the planform and the diaphragm forms 


Se 


tials) may chen be determined as for any other purely supersonic 


surface, 


The illustrative examples shown below are carried out with 
the PiC-method. ‘The ViC-method follows an identical procedure 
with these exceptions: 


(1) 


(2) 


(3) 


- we ee ee 


One aust employ the VIC's ($, Ron te hod instead 
a si 4 


of the PIC's (Cy = CPP ae aa) Both of these 


7 


quantities are defined in Appendix A. 


For diaphragm regions near ¢ subsonic leading edge or 
a side edge parallel to the flow, one must impose the 
boundary condition of 4¥=0O at the centers of the 
boxes rather than AB=o. 


For the Ciaphragm behind a subsonic trailing edge, the 


true boundary condition of AP =0 dictaces the use of 
the relation ( x, ¥ being behind the trailing edge) 


1g fX~ XR.) 
Aybiy)= s9¢(ranye BE? Eq. (4.1) . 
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this expression is derived from the os relating 
the pressure to the velocity potentia 


Ap = gif hape ¢ 60) } Eq. (4.2) 


and the condition that Apuo behind the treiling edge, 
on the Dia 


IV.1 Determination of the Downwash Distribution hra 


eg on 


If the downwash over each disphregm box is assumed constant 
the diaphragm region should be represented by those compiete boxes 
of the grid system having their centers ineide tnis region, Im- 
pastes the appropriate boundery condition at the centers of these 

oxes, the downwash distribution can be determined by “sequential 
solution" of the resultent equations. As an example, consider 
the rectanguler wing shown in Fig. IV.1. For this case, the dia- 
phragm boxes are (22), (3.2 ), (42) and (37), 


ey," 


Fig. IZ.! (tlustrative Example for the Treatment of the Side Edge 


One has the following equations; 
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og eee 


=, 


Eqs. (4.3a-d) 


it is noted that in each equativn, there is only one unknown, 
provided all the eco | equations have been solved. This 
ailows a simple Sequential calculation and avoide the necessity 
of solving a sim ltaneous System. If a matrix formslation is 
used a8 in Ref. 8, there follows 


ii, 
4% 
Ge EO 8S) Co aE Bs 
| - i, 4 
Gj Ba O_O 6 6 Sy atel, . ad iy : 
uy 
C,, io G. 0 0 o : Cy 0 es “ re (ed Oo 
Co. “al ) 
Cy Cy G2 G Ge Gs! Gel, Ge G is, 
Wa 
Eq. (4.4) 
or equivalently, - 
Ge? 0 9] Gy 
Cy Go 9 9] Ge 
40 ’ Coo a he 
0 oy Gy G, fm G Eq. (4.5) 


WADC TR 56-97, Part 2 


CONFIDENTIAL 


It 1a known frem steady-state solutions that the downwagh 
at the side edge, when a proached from the diephragm Eageem, 
exhibite « square-root singularity end vanishes at the oremost 
Mach line. The assumption of constant downwash over each box 
of the diaphragm therefore introduces errors in the pressures 
over planform boxes near the side edge. If e Large numbar of 
boxes is used, these discrepancies are confined Co & smell region 
near the tip, and the resultant accuracy for the tote) generali- 
zed forces ig satisfactory. The process just described ie the 
one Suggested by Alternetive 1, Section 131.2. 


To improve the airload estimation, es eclally when a small 
number of chordwise boxes is taken along the side edge, the ef- 
fect of the Singularity can be accounted for more accurately in 
the following manner. 


Referring to Fig. I¥.2, in Cartesian coordinates, the re- 
ce point is (%% )while € > are running variables rep- 
resenting the influencing point. ‘A characteristic coordinate 
System rs with origin at (9) is definad such that the axes 

rand Ss are, respectively, the right and left forward Mach 
lines emanating from (“9 - The new coordinates * and ¢ are 
related to x,y £ and 7 hrough che equations 


r= el (x-8)-Oy-2)] 
$e Bi! (x-E)+ BGy-9)] Eqs. (4.6a-b) 


If the receiving point is on the aide edge, 1.@., J~O , one may 
assume a downwash amplitude distribution 


G (iis) 2 ear Eq. (4.7) 


Vd Vrs 


which exhibits the Proper behavior at the side edge (r*5 or 
g=0 ) and vanishes on the foremost Mach line rer. Here d 
ie a typical length and #7; is the "strength" of the downvash. 
For convenience, a2 will be taken as half the diagonal Length 
of a Mach box. 


It is more convenient to locate the origin of the coordinate 
System at the point for which the pressure or potential differ- 


ence if to be determined. By coordinate transformation one ob- 
taina (see Figs. IV.3a-b) 
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ALONG THIS SIDE EDGE ras 
ns erin arian ¥.7) 


(r#0,680) 2. 


x,€ 
Fig. IZ.2 Characteristic Coordinate System 


4“ (4,8) = ee 


Vr+g,~$ 
hq. (4.8) 
when the origin is on the diaphragm region, and 
“Fer 3) sa wi, ee he Eq. (4.9) 
0 Md pags 


when the origin is on the Planform, 


It may be shown that for steady flow the pressure coeffic- 
lents associated with these two cases are 
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(b) 


4 ~~ 


Gr Lownwash Disiribution Near a Side Edge 
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(il =f =O" fr fy use 
om /4- + [F=f Fig. West 


Eqs. (4.10a-b) 
Similarl 


y. for the velocity potential coefficients, one hes 
(4% being the streamwise dimension of a box) 


PEE $e a he ms 
é 


SG GIB)" 6(G)"s af fe Fy a8 


Eqs. (4.10e-d) 


This refiuement may be introduced in the Mach box scheme. 
To Lllustrate its application, consider the Simple problem of « 
rectangular wing in Steady flow at an angle of attack & . 
The downwash at ail points on the planform is //~ . Referring 
to Fig. IV.4, assume the following downwash diatributions; 


(1) Singuler distribution of strength «<{ over the region 


Cac, 
(2) Additional singular dietributions of Strength «:, , 


wr Por over the regions ¢ ek ehé 
respectively, nae: 7 ae 


’ 


(3) Constant downwashes of 


Strength «5, and «;, over 
poxes (5/ ) and (4% ), 
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Fig. Iz,4 lilustrative Example for the Treatment of the Side Edge 
Including the Effect of Downw ash Singularity 


Satisfying the conditions of sero pressure at the centers of 
diaphragm boxes, one has 


fy27 9+ CG, (te) Co (fen 1) tty “0 


P= Oi Cae (Ue) 2G ni) G, wz, = 0 


As “OO; Coe (Uk )+ Cy, (Ou) + Cy, (Ue) + CG - 2) 4, 
* Q'G = 1) ti, = O 
’ 
Eq. (4.11) 


As in the previous example (Eqs. (4.5)), « Sequential solution 
or a matrix solution of the above set of equations yields the 
desired downwashes on the diaphragm. 
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IV.2 Determination of the Pressure or Veloeit Potential Distr;- 
butions on the Planform a snnt_Locentna) Distrt 
eA eet RE Rone romtNanerA aeetoomeen 


Once the downwash distribution in the diaphragm region is 
determined raps del, the previous section, the pressure or 
the velocity potential differences at the control points of the 
pranfore are next obtained. If the downwaeh over each diaphragm 

ox i8 assumed conatant, the procedure is identical with p- 
pendix A [Eq. (A.6) or Eq. (A.8)1, where only the preseure co- 
efficients Cy xz or the potential coefficients Ps x aséociated 


with unit constant downwashes are involved. When the singelerity 
in downwash is to be accounted for, the method ig similar but 


additional tabulations for the coefficients Ang GC” (ord 


€,” and es ere necessary, For instance, the pressure at 
center of box (44) ia (Fig. Iv.4) 


fay | he lob, * Coa Goa? Gat Cae # Lag? Ce, 
+ Cart Ci, Cz, + C3, + G) 
is GPa 5 BoD) Mee GB nt Saa)] 
Eq. (4.12) 


IV.3 Aerodynamic Influence Coefficients Associsted with Sin uler 
wives stribut ions ne 


The general expressions for the coefficients ae : nd 
@ and db,” are as follows: 


OP bis Se | Lea- Le a} “Mg La (28) 


ih 
Sead Batis. 
(3 


- 


G'= bie E{le2)- Peal} -2 [ER (a) - iat 
6, 


Eqs. (4.13a-b) 
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@ ~-4(B)" {Beor) - Be aa re ez) 


Eqs. (4.13e-d) 


where 
w is the cizcular frequency; 
a is @ typical length feat to be half the diagonal 
length of the Mach box 
a4, % are lengths along the characteristic uxes aa shown 
in Fig. 1V.3; 
G st Le (wo a 4 (H4 
a6 se Mild 
& 
tJ 
% i 
= % 
or = e 
and 


*. i. n ati wl ‘+ 
} my) 
L,(69) er [ ipa ae EF Vrs ad Sosy vrs) dsdr 


7 cu ma Mr+8) “0 
£, (68)= “mat | [mee ae oF “tet Vrs )dsdy 


Eqs. (4, l4a-p) 
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For the steady Case, Eqs. (4.13a-d) reduce to the simple 
forma of rg (4.10a-d). In general, however, the integrations 
for J, and - must be carried ont numerically, When Comput ing 


f”, one need not tabulate J, and J (or fa and / ) inde- 
pendently ae in the case of aaa By suitable transformations 
of varieablea, the double integrals reduce to the single integrals 


L, (ga) me ae x Mees melee) Ae [evc]) 
+ Zeng ( prelate) J (ele +) 
iedd 


HZ al eB) T (¢fa%el)} az 


fener ' / ‘ 2 
Whoa) 0 | oat Se se xsi re peta 
Ve og 
+2; poe ( Sea eg SZ, (fe*a] ) 


ico" alte, kor) 


£qs. (4.15a-b) 
where J, is the Bessel function of the first kind of order 4. 
The infinite series in the integrands are rapidly convergent for 
sw terms need be retained. it is worth noting that the 
integrands are functions of »* » O is an integer or a fraction 
of integers and = ig a fraction of integers and less than / 


The numerical method Suggested here is a modification of 
Gauss’ quadrature. It can be proved chat the incegration forma 
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[Sede = 3314 467) ba. (4.16) 


is exact if //2*) is a Polynomial of degree ( 27-/ ) or less in 
eo provided the a are chosen properly. H,’s are the 


weighting factors asscciated with this set of 2; (eee Ref, 1). 
For V=S , one obtains 


Z* = 0.022 (63, $67 H, = 0.295, 524 216 
z/ = 0, /87 B31, 57% bh = 0.269 266,739 
Zy = 0.461, $97, 34y Hg = 0.2/9 086, 348 
2," = 0. WE, 334 658 Hy = ONY YE) 36/ 

Z, = 0. 948 483 9/0 He. = 0, 066,671, 338 


Eqs. (4.17) 


For low to moderate frequencies, the integrand of Z, can be 
satisfactorily represented by the polynomial 


4 ‘8 
F(z*) = Ata, 2 ng, aig foes +a, 2’ 


which passes through pointe z,, ... ‘Ee 5 and the integration 
for J, accerding to Eq. (4.16) will be adequate. 


For £ » the situation is somewhat different since the 
range of integration veries with 3 -_ However, it has been 
found by some trial calculations that Eq. (4.165 can stili be 
used with the following modiffieation: 

ae vA Ue a 
Lf, a et [ gta 7 Btn {Qle*)} de Eq. (4.28) 
ver. 5 ee - os © Ue 
" ef UL ze Oe ide - y5°tt ee Fey] 


where G(z*) is defined to be 
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Mle)- {4 I Cpe 2 vite) L(efe*a]) 

+ £ COL eee lefeee)) 

NEL PME) el) wher ose 
= 3 when z°<F 


ten Eq. (4.19) 
By numerical integration £ becomes: 


‘geal 4 a 
F) xe? >) fm meg Og Me ere 
LbF) = @° Lb lge “alle;)~ Jr we ‘} Eq. (4.20) 
For higher frequencies, the integrand becomes sinuous and 


caseer degree Polynomial must be Considered, Reference | dig- 
8 situation, and outlines the Steps in devisin a 
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APPENDIX A 
PRESSURE AND VELOCITY POTENTIAL CALCULATIONS 


A.l Pressure Influence Coefficients for Square Boxes 
Senne recs Lor Square Boxes 


Although che present report concentratee on applications of 
the Mach box scheme, influence coefficient formulas are listed 
first for the Square boxes, since extensive tables of these co- 
efficients already exiet. ‘Tha preesure difference Sri at the 


center of the receiving box (2,7) due to a constant downwash 


We” at the sending box (12) is (Fig, A.1). 
us ‘ ad >i ae é 
= anal ty teary Li See 
Ro, aia Sidi Wea Cig, + Ly g) Eq. (A.1) 


where Vamn-y , Gem-m and Keil 4, the associated 
poets influence coefficient, As depicted in the figure, the 
ntegers and ™ indicate that to reach the receiving box one 
moves 7 boxes downstream and ™ boxes to the right from some 
arbitrarily designated reference box (2,0 ), rr 7 and 
represent the number of steps downstream and to the right, 

respectively, that one moves in going from the sending box to 
the receiving box. In supersonic flow 227 always, 


The total pressure at the center of (27) is therefore 


a Yo 7 -~ er 
= = ae - (A.2 
fae thin OE EH Regt Ty,) tw ta 


where the sumeation extends over all boxes ( uw ) that lie par- 
tially or totally inside the forward Mach cone emanating from 
the center of box (rm). These coefficients ame tabulated in 
Ref. 2. Only positive values of /* need be considered since 
by symmetry, 


tit = Ko tea Log Eq. (A.3) 


A slight error appears in Eq. (7a) of Ref, 2, and ie cerried 
through in the tables. The correct expression should read 
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= a2 m4 
Re £[-1+ 22 (fy- Hee fy] Eq. (A.4) 


Furthermore, the accuracy of Ke in the moet extreme cases 


(f.e., large 4, at low Mach numbers) is no better than five 
units in the third decimal place. For these high frequencies 
and low Macon nuwbers, thie is not serious, since it is known 
thet the accuractes for the rows Y>2 are much poorer. 


A.2 Aerodynamic Influence Coefficients for Mach Boxes 


For Mach boxes, the pressure difference at the center of 
the receiving point (7) due to a conetant downwash ee 
at the sending box (2%) is defined as follows (Fig. A.2): 


wee PY ett ; 

SO emennataeeetetrion Pr E e A.5 

aa ae Gy, (%, +4 — q. (A.5) 
This definition is more convenient for tabulations, and for the 
steady case R, z 418 independent of Mach number. Accordingly, 


fe 
the total pressure at the center of (4m) ta 
ue CT wt 
mL) Pm PLY tut sy — Eq. (A.6) 
ho yy Lp (+x) 
The velocity potential difference at (%™) due to a con- 
eteant downwash %, evr at che sending box is 
"AO 6 Me 1G? gM 
Ay” oe (MW ee Eq. (A.7) 


The total velocity potential difference at (7,7) is 


tat i) ) 
AY = yy" 6, — Eq. (A.8) 
es ie han ge OG Org Sg 


Again by symmetry, 
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Re a * iS 


= 
? “0 fb 


» Oe te J 


Ei st 

DH tp A | Ftp 

1, aa “* ae = Ws = “ae 

Eqs. (A.9a-b) 


@. Exact Expressions 


The exact expressions for the PIC's and the VIC's are now 
given. These were derived by using the series expansion for 


the velocity potential given by Watkins of NACA (no formal refe- 
rence available, see Ref, 1). 


“me a 
Rrit,~-e TR) £ELME), — sizec 
Fegr iff, m4 be OU bani t DE] tog onl aoa) 
; re . i 
ERY © Aedes 3) St ann SE 
B-¥ ray re 


Dm fe > 


to he lasts 21 Goes) bron ss.)] 


-hiy- } : 4 7 n. 
+ £6 “her sie, + E05, Ge sonfrsints 


2o-/ ray J 
; 2. a e; : ; Pie 
- 04 oN insy + SCI CE ambrsin Lyk, 
*t V>r1 
no 
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® 5 f- Poi H(3 ri fy ee Pe ieee) 


Aor ayrl 


fa] 


: Py ATE opin _ sin" Sit!) 


ands 


LE te a) snes e*)- sulersin” 4] 
+ hey ye A yin we. FY a i pores Fi anor gH 


oh! 


ty amg 


Eqs. (A.1l0a~-d) 


Aref” » -fin¥) , vam ee 


aft 2 f?- of? BX (Lachspow'eet F Fr, ye uy r Eijsmbersn” 2 1 evs zi) dy, 


Vm >/ 

Revs. y ray 2 AL y)s0't os BOT then lrsin LAX, 

Vr, Aro 
OF lx Ry oi 1ae-l ot 2a 
ari 8] € tan eon se) 
ar 3 
Bl] Gryfon len’ Be) selene Be 

U>h>/ Eqe. (A. 1Lla-d) 
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where 4, is a reduced frequency based on the Streamwise dimen- 
sion of a box 


Cae 
ee Eq. (A.12) 

and J is the Bassel function of the first kind of order 2, 
The function f 1a cowmonly encountered in linearized unsteady 


Supersonic flow (cf. Refs. 9, 15). Although the expressions 

Seem rather lengthy and comp tcaced, the task of ta ulating these 
coefficients is not difficult, because of the rapid convergence 
of the infinite series. (Ic is assumed thet hig “speed com- 
phen will be available, The sine functions in the series ma 

e written in terms of the Chebyshev polynomials S, (x) (Ref. 10 


Sinfarsin™'x) = (1) Jr ae 2X) Eq. (A.13) 


b. Approximate Expressions 


For Mach boxes, the approximate expressions for the pres- 
sure influence coefficients are next given. These involve the 
same degree of approximation ag those of Ref. 2 for square boxes, 


) 8 ti hy (orn £3 : 
Sled EY AE CEUB) , 2-0 
“oat Lol [Epon nfo ah Poon yh) 
vil MOG ontop aE, ponte pi 


Vim jt 22 


Gre g f } AenfhloloaIG - Lost ope coal, fivg = di ®) A 
tebe Mesh RM, vit Ecol os AG, 
rb oti eM tlm le TEM} 


U>f 22 


Eqs. (A.l4a-c) 
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where 
- =) * 
Ie =~ 20d bn. 9, = IE ees £ =H» 
Ba ads) ‘¢ 28 f o 
i : et Po2, =o 
Eh, cosh Level) - cash (ap-1)+(0+ he - (OHH, , 
Sal = 
- wt gem ? x 
Mop on" $F, Ty 20, Tujeré 


G5 = -(o- fjeoh” ZH (wig te, , I 


r “= o ~ mn Sad 
ot Bae) ov Bsae! v 2t=/ ae ad 
Mw cas SEH! eng ait A mcs 82! Atl . of 
GA FBS B61? pe B55 OS BS My 
Vee? 


- ley. “ev-1] fe. + 2TH ov gh 
Cp fir $)fush ot ef é fous oi-i pide fa! Pos>/ 
+ (Deg) Aer ~ (0-4) ty, 


ise interpretations of as , cash~ 
2.l0a-h), Ref. 1, have been applied.) 
For the boxes in row @=/ , a subdivision da sesemet 
which is similar to the method of Ref. 2, Reforring to Fig. 
A.2, each such box is divided by a finer grid of sxs elements, 


» etc. given by Eqs. 


‘ Pot z . ' 
GL 2 Oarsl, rt Diep L ji > vo}, fewo 


Gyre = 22? (Cpr Fen dy 5 Fae eA 


Eqs. (A. l4d-e) 


where the subseript indicates that an rity Zz depend on a re- 
duced frequency argument equal to one fifth of 4 and are ob- 
tainable from Eqs. (A.l4e-c). 
Daa — 
As may be seen from the definitions, Aya ; ys and 
Bie are independent of Mach number, in contrast to the 
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corresponding coefficteate for the e uare box scheme. The Mech 
box tabulations are therefore such simpler. The values of 

a ; He and 6. are given in Table A.l for = 
from 2 to 25 and for Z from o te 25 . Note that onl 
values of 4< 7 nesd be considered. For the ¢ eady case, the 
above formulas are exact. A short table of va (4-0) ie 
also included (Table A.2). ' 


For the velocity-potential-difference influence coeffic- 
ients, the following approximate formulas apply: 


te ; a » a it) = é 2 i ta 
Moot le ef Sh) Sacer ye inlyel). ets RY), oom 0 


yy rt" ; or = 
Rovefe FW 2 (cosvh, -<s 98 )cos (Be la=®) Gl, 7 I Ves P2 


oy 


YM tA Lx 7) w) 4 21209) ree \ = 
ror J, Ailes ht 8 BE Ons LO) ee 


B23 3 , Basems ( 


Sie y a a tar, 
Cae he ae ptr a CAM Me. rig) 
o FF 


Eqs. (A.15a-d) 
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APPENDIX B 
DETERMINATION OF DOWNWASH DISTRIBUTION 


In the determination of the generalized forces for use in 
the assumed-mode “pproach of flutter analysis, one wust first 
obtain sither the pressure or the velocity oOtential distriby- 
Clon over the wing due to the downwash which is associated with 
@ach of the assuued free vibration modes. The downwash ats, yt) 
Forresponding to an assumed wode Af,y) in simple harmonic motion 

complex representation) may be expressed by 


‘ae ewe 
als yt) = Ot fla pe é Ley} . Eq. (B.1) 


The downwesh is a function of both the deflection and its stream- 
wise slope. If the Sraieetity of each such mode with respect 
to the streamvise coordinate x ig known analytically, then no 
difficulty arises in determining the downwash at any desired 
Point. For example, in the case of « beam-rod type structure, 
where there are no Chordwise elastic deformationa, £¢ #) ie at 
most linear in the variable x , However, for a piece Cype 
structure, there exist chordwise elastic 

muse resort to graphical differentiation or construct an analy- 
tical approximation for Siy) in the x-direction at a given 
spanwise station in order to obtain streanwise tlopes. 


In practice, the flutter mode is azeumed to be some Linear 
complex Combination of the first few free vibration modes, For 
a pla 


iecrete points are usually obtainable for these modes. The 
main task in evaiuatirz the downwashes is to find the stream- 


the deflections at these given points. If, at each spanvise 
station corresponding to the Centers of each column of boxes, 
there exists a sufficient number of pointe a+ ach deflections 
can ba computed, then Simple graphical rocedures cen be re- 
forted to to obtain the downwashes at the centers of thea boxes, 


‘Otherwise, the usual procedure is to find suitable polynomial 


expressions for the deflections, Then the Slopes are obteinsble 
by proper differentiations of these Polynomials, 


B.1 General Procedure for Finding 4 roximate Anal tical Expres- 
sione for the Detlections A = 
CN naa ee CED 


Consider a cantilever rectangular wing on which the deflec- 
tions Oy) are given at min points 43 shown in Fig. 8.1, 
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Given the defiections at | de Hirt’ Z% at any chord- 
wise location % one can find a polynomial in for the deflec- 


tion Fiteoy) . For thig cantilever case, {f the polynomial ig 
given ly 


IO)" CaP tl gie + hag a ayy, Eq. (3.2) 


then one has 


Bl)" Ge > oe Li rhyme 
! 


q 
#6 be) bates hase + batt 
Eqe. (B.3) 


g(x,y) 


nm yj Yn 
Fig. B.! Typical Sponwise Deflection Curve 


Solutions of the above eet of Simultaneous equations yield. the 
constants »* é, si 4 Gy, + Analytic expressions for other 


chordwise lecations way be obtained in a cimilar manner. In 


Seneral, the deflection Over the whole surface way thea be 
written as: 
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FO, y) = Abele) + bl Ip ly)e oot Sv) hl 9) Eq. (B.4) 


where the functions 4, (1) are Polynemials in 2% which satisfy 
the relations 

4; (x¢) m/s» Bi (hp) =O when hy? Eqs. (B.5) 
Again one may find polynomial expressions for Aix) » The de- 
gree of £(x) will depend on the number of chordwise oints. Once 
the general form of the deflection [Eq. (B.4)] is o tained, one 
can easily determine the downwash at any desired point according 
to Eq. (B.1). It should be noted that either a graphical method 
or a procedure similar to the above for finding the downwash is 
necessary regardless of the type of grid system used. 


Tt 18 worth mentioning at this point that when one finds a 
polynomial which passes through a given set of Points X,,5) hy 
(or yj x ), one can be sure of its reliability only in the 
range 4) 4X8k_ (or vey, ), Beyond these limits the poly- 
nomial may deviate considerably from the true picture, It is 


larly, the centers of the boxes nearest the tip should not be 
much beyond y, . To insure this in instances where the deflec- 
tions at points near the edges are not given from normal mode 
data, some artificial means must be devised. One such means is 
ee ri eee of the deflection along constant x- and con- 
rion Segal So that a reasonably smooth deflection surface 
results, 


B.2 The Triangular Wing 


In essence, one might consider the above as & useful scheme 
for finding deflection polynomials for other planforms, such ae 
a triangular wing. As an example, consider the oversimplified 
case of Fig. B.2; deflections are given at only six oints, On 


t is necessary to have 
artificial, reasonable estimates of the deflections at points 
@,&,¢ which may or may not lie on the plentorm. However, such 
calculations demand cut and trial cross-plotting procedures of 
deflection vs.x% for constant values ofy and of deflection vs. # 

for constant values of x with arbitrary extensions of curves 
whenever the data are lacking. The following procedure (which 
can _ also applied to the rectanguler ving as well) is sug- 
gested: 
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Fig. 8.2 Itlustrative Example for Finding the Analytical Expression 
of the Deflection on «a Cantilever Delia Wing 


(1) From the aforementioned cross-plots, find at the given 
sets of x, the deflections at the spanwise stations correspond- 
ing to the centers of various columns of boxes. When changing 
Mach number, the spanwise location of the boxes will not change 
if the span-dimension of a box 1a kept fixed (as suggested in 
Section III), 

rs Knowing the deflecticis at a discrete set of points for 
each desired alte station, one can proceed analytically to 
find polynomials in x for each such spanwise station, 


(3) Using these polynomials one can differentiate to find 
streamwise slopes. 


In the above steps, one may se St any existing additional re- 
Dg Peat on the polynomial which would supplement the data of 
eflections, For inetance, one can utilize the boundary condi- 
tions for the cantilever wing 


Zed F ge (7) jar af x's al y=9 ’ Eqs. (B.8-9) 
Another condition is that the trailing edge is free, so that 

2 z x Fp age oP ee 

(2. Pr vgr)e = at traili ig edge, Eq. (5.10) 


Ag it stands condition (B.10) is difficult to use, so that one 
may take this condition approximately by setting 


3 2 O at trailing edge. Eq. (B.11) 


This is justified, since for built-up beams the effective 
Poisson's ratio vw is small. 
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B.3 Lagrangian Interpolation Technique 


In cases where, for eech mode, both the deflections and 
their slopes are given for a set of discrete peints, it is pos- 
sible to find the deflection and ite slope et any other point 
by direct interpolation. Then the determination of en analytic 
expression for the deflection ig no longer required. 


Given a function { at (24/) equally spaced points, it 
is possible to pass through these points @ po ynomial of order 2 


fih= Kd) Peey)t- +h NL mo)e +k A)SE,) Eq. (B.i0) 


where 2= (x-%0.)/A and 4 is the spacing between points. it 
may be shown (Ref. 11) that, for instance, 


Kam -$la-a*'] , Kaa s-a*, Kd= $/aea¥ ; for N= | 


Kiln A*-247-A 424 


- AA bat 4) 
ue K.,(A)~ ede ba td 


? 
MAA HAY 4A 


#542 = 
K(j=»Aasa ee Ki, (A) . 


AY, 24° A 2) = 
Ky ta)= -- aaa j fac N 2 ' 
Eqs. (B.11) 


Tables are available for the “s vg. A_, but if A is irra- 
tional it might be easier to use these formulas directly. In 
cases where the given points are not equally spaced, one still 
can devise epprep: iete interpolation feeeatar es which ere simi- ; 
lar but somewhat more tedious. To il uétrate the application 
of the above, consider a spanwise etation of the wing having 
eix equally spaced chordwise stations 1, ... , 6 where the 
slopes 6; are known (see rie B.3). It is required to find 
the slopes at other points pecan as shown in the figure. 
Assume that the slope at point 7 may be adequately computed 
using a three-point interpolation, i.e., jms. 


Sat 
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Fig. 8.3 Logrengiun interpolation for a Function of One Variable 


Then one can write for the deflections at the new points Li F 
6 = KA)8 + Ke (Ar), * Ki ly) Os 
! 


a = K., Ay) G+ *(46) 6, - K, (4g) Gg See. (B.12) 


or in matrix nNotaticn, 


{o,} r Leaf; Eq. (B.13) 


where {6r} and { 6, f are column matrices, 


5 a oh he 7 a o ° ] 
0 K, fg) K, (Ay) Kia) nd e 
fica iJ ra) 0 Kj(As)  KelAs) = Ks) O 
ee 9) 0 ro) K.(Ag) Ke) K, (ar) 
re) 0 o K.,(Az) Ke (Az) K, (ay) 
Eq. (B.14) 
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and 
_  ABn-Xe “EA re rs ee 
Ay oe ag Tilt - A, 4 ty AL E18, (8.15) 


In @ Mach box grid system, if the Mach number is changed 
while caeh ing the spanwise dimension of each box fixed as 
recommended in Section I11, the control ointe will change posi- 
tion in the chordwise direction only. Therefore, the type of 
interpolation matrix illustrated above is sufficient to deter- 
mine the slopes ac the new control points in terms of the known 
slopes of the original set of fixed points (for the initial Mach 
number). In instances when there are fewer than three points 
for interpolation, such as near the tip of a delta wing, a 
lower-order interpolation must be resorted to. 


BM Lagrangian Interpolation for a Function of Two Variables 


To obtain the proper interpolation for a function of two 
variables, the above technique may be used repeatedly. ‘two 
such examples are: (1) when the force-displacement. structural 
influence coefficients C,, (which is dependent on the two con- 
trol points </ ) are td be found at different chordwise points 
for given spanwise stations (see rest i and (2) when che deflec- 
tions and slopes are needed at points %,¢ ) for which x ¥% do 
not lie on the spanwise or chordwise locations of the discrete 
points with known deflections and slopes. (In the previous sec- 
tion, the ¥ -coordinate of the new points coincided with the 

-coordinate of the old points, so that the interpolation 
was in the stream direction only. 


As an illustration of case (1) above, let <(or ¢ ) define 
the original set of points and Z (or / ) define the new get. 
Bearing in wind that ¢ and< (or ¢ and / ) are restricted to 
the sawe spanwise location for this example (Fig. B.4), one ob- 


tains 
[ CyJe [KA] [Cy] Eq. (B.16) 


and 


lo, ) = LG IKuy]” 
Kaa 


Eq. (B.17) 
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where { J” indicates the 


transpose of the matrix, Re- 
ferring to Fig, B.4, 


STATION ® , 
STATION @ 


Fig. B.4 Lagrangian Interpolation for g Function of Two Varicbies 


K., (A,) K, (d;) k, (A;) 


Lea} K.,(a3) K(Az) fe) 
Eq. (8,18) 
[6 Ge Ge] [Kay Kt) 
" 
Ic}- Cy Cos Cx ) [K(4)] = K, (Ag) K,(As) 
“j 
: Cyy Css Css kK \ é Ay) K, (Az) 
Eqs. (B.19a-b) 
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Therefore one obtains for the new set of influence coefficients 


Ic, J 


ae )  K,,fA3) 
KO) Kay Kl) = ta 


‘ GC. Ge Gel {KQs) Ky) 
KA) Ky) Kia] | 7% 
Cu Ki (;) K, Ge) 


> 
> 


Eq. (B.20) 


The above procedure may be repeated for all other epanwise 
stations © and @ to obtain the complete set of influence 
coefficients. 


For case (2) above, the situation is somewhat different. 
It entails interpolation in one direction, say ¢ -direction, 


rise * constant, followed by interpoieation if the other 
direction ( x -direction), keeping f constant. 
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APPENDIX ¢ 
RECOMMENDED INTEGRATION TECHNIQUES FOR GENERALIZED FORCES 


Once the downwash distributions are determined ag Sigs | to 
Appendix 3, the racommended box schema may be used to find the 
pressure (or the velocity ge tld distribution aseociated with 
each assumed mode. Then the generalized forces are obtained by 
integrations over the planform of the pressure (or the velocity 
potential) distribution weighted with the various mode shapes. 


The pressure distribution along a streamwise 6trip in che 
supersonic region of a planform i¢ continuous end smoothly vary- 
ing. To evaluate the c ordwise portion of the double integral 
representing a gencralized force, a simple integration echeme 
such as the rectangular rule will suffice. However, whan a 
6treamwvise strip includes a Segment of the mixed wingtip region, 
the chordwise peesture distribution exhibits a sharp drop end a 
n slope across the Mach line. If many boxea 

were included along this ‘strip, the segment between the Mach 
line and the succeeding downstream point (i.e., the segment con- 
taining the sharp drop in presence distribution) cortributes 

the integral secross the entire chord. 
Hence the error introduced will be small if the rectangular rule 
le used. However, if there are only a few boxes present (say 4 
six or less) across e streamwise strip in the mixed re ion, the , 
contribution to the total force dus to the sharp tg nm pres- 
sure distribution may be appreciab:s. A refinement in the chord- 
wise integration procedure is then recommended. 


C.1 Refined Chordwise Integration Scheme ( PIC-Method) 


Referring to the planform of Fig. C.la, a typical chordwise 
pressure distribution for section A-A i8 shown in Fig. C.1b, 


From points ¢ to 2 and 3to d » where no sharp drop oc- 
curs, the simple rectangular rule may be applied. Between points 
Zand 3 , the following integration technique is suggeeted: 
let the pressure distribution between points 2 and 3 be 
represented by: 


Aix) = (1- / ) +k (V5) Eq. (C.1) 


where ¢, is the length from point 2 and 3 . This expreasion 
yields the correct pressures at 2 and 3 , and has an infinite 
slope at 2 . The infinite slope condition is imposed because 
it represents the correct behav or, as may be derived from 
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Sa 


LEADING EDGE 


TRAILING EDGE 3 4 «68d 


(a) K (b) 
Fig.C.l Typical Chordwise Pressure Distribution Near a Side Edge 
(Steady - Stats Condition) 


steady-state results. srrnpresing A(x) vecween points 2 and ’ 
3, the formula for the lift in that segment is 


af ~ LO" ply 
= ¥ (href) 


The above integration in conjunction with the a eae rule 
for sagwents other than a-3 , ylelds the Lift per unit span at 
A-A 


ge 


Eq. (C.2) 


a¢ 
f- / Plx)dx 
“Gh+ G(B)+ Bina) B)r bbe bp 
= G(b+ ZA+Zhrs +A) Bq. (€.3) 


This refined integration formla is very siaiisr to the rec- 


eapgucer rule, and its use is no more difficule than the simpler 
rule, 
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In a similar wanner, if the moment about any axis %, ig 
represented by 


x, 7 foe 
Om | OPP) AE Y-m dean. (e.m 
Xp , 
there results the moment expression for the segment 2-3 


dm Fb Blk re shen, bn)ef (bdar Bt, $x.)/ Eq. (C.5) 


The above procedure cen be extended for the calculation of 
the generalized forces. If the weighted integration between 
points 2 and 3 ig given by (assuming Linear variation of the 
weighting factor between 2 and 3 ) 


= r fia ae PER) Lf [RR PN fy) 3d 
APL PET Bl) AE ROE SOEI wa. (0.6) 
2 % 


then thie approximation yields 


Ar = 6p (afte) bfte))+ h(a Male $/t))] Eq. (c.7) 


Ie should be emphasized here that thie refinement must be used 
only when there are no severe fluctuations of the pressure 
across the chord. This fact iimite its use to Planforms with 
supersonic leading edges. For planforms with subsonic leading 
edges, the whole planform is influenced by the diephra region, 
and then one must take a sufficiently large number of oxes and 
use the wectangular rule. 


C.2 The Determination of Generalized Forces PIC-Method 
eee’ Forces (PIC-Method) 


A generalized force is given as a double integral over the 
planform of the pressure of mode < weighted with the deflection 
function of @ mode 7: For instance, in dimensional coordinates 


xg 
Ay, * LY A 06, yd dy Eq. (¢.8) 
s 


The integral over the chord at any #panwise station bo is 
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= AB Hm, aT (Am, 
Jo priftendde = Sh a564) 4h bow pf ( w} ree 
where 4, (4) sra the adjusted integration coefficients of the 


prcvious section. In cases where the rectangular rule is suf- 
ficient, the 2, (%) Will be equal to the box dimension / ex- 


cept for boxes at the trailing edge where some sort of an area 
correction is recommended (cf. Rule 4, Section III.1). Im addi- 
Petal @ spanwise numerical integration yields the generalized 
orce 


NM 
Q, SZ) 4h lte) {AB (rm 4) 4 Crm by) Eq. (C.10) 


where d,'s are the epanwize integration constants. If here the 


Tectangular rule is used, the 4,5 are all constants and are 
equal to the spanwise width of the strip (4 /e). 


C.3 The Determination of Generalized Forces (VIC-Method) 


Using the pressure-velocity potential relation 
Me: oa 


the expression for the generalized force (Eq. (C.8)] may be put 
in the form 


q = sd [fd ARO Mila ylbucy + sey, (fo CLE i yib Eq. (C.12) 
S Ss 


A partial integration with respect to x for the second area- 
integral yields 


Uh | 
aad FL agiug [bf tay é Chow dx dy 


+0 f CAC DIAC Eq. (C.13) 
Shen 
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where the subscrip: 7 refers to the trailing edge at & + Here 
one may use the rectangular rule for both the x- and g~ inte- 
erations, with the result 


§ lf lH nif fit Liaoetcayli fla) He) ] i Eq. (C.14) 


Here 4,=4,/e and Z,(y,)=# except near the trailing edge 


where they aust be adjusted to comply with the area correction 
rule (cf. Rule“, Section III.1). 
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APPENDIX D 
FORMULATION OF THE FLUTTER PROBLEM 


It has been common Practice in the past to Qppt. bp a« flute- 
ter analysis as a Series of related but separate problems. The 
flutter mode ig assumed to be some linear combination of the 
firet few free vibration modes of the structure, These vibra- 
tion modas are calculated by whatever method ig feasible. Next 
the aerodynamic Problem iz solved to provide the neceesary 
generalized forces associated with the aseumed modes. If two- 

imensional forces are to be rrerees aa in cases of large- 
aspect -ratio surfaces (Ref. 12), this step is considerably sim- 
plified. With these results available, the flutter equations 
of motion may be formulated and solved in any of a number of 
possible ways. This Procedure has been very successful for plan- 
forms which are rigid in the streanwise direction, l.e., the 
streamwise slope at any epaneiae station is constant for any one 
of Che assumed modes. A bean-rod type structure falls in this 
category. Exparience hes ahown that even rather approximate as- 
sumptions as to the mode shapes and knowledge of the free un- 
coupled vibration frequencies of these modes usually yield satis- 
factory solutions, However, when Plate type structures, such ag 
deite wings, are to be analyzed, accurate representations of the 
assumed modes become necessary, 


Another approach is a direct incegral equation formy lation 
for the flutter problem in terms of aerodynamic and structural 
1ufluence soefficients (Ref. 13), but this method does not offer 
all of the advantages in simplicity obtained from the assumed 


D.1 Equations for Bending -Torsion-Aileron Flutter 
lf ¢ 


he structure to be analyzed is of the conventional beam- 
rod type, the formulation pe to the present report is very 
Similar to the method presented n Ref. 12 with the sole ex- 
ception that the Seneralized forces are evaluated by the aero- 
dynemic influence coefficient method (the Notations are those of 
Ref. 12, except that the spanwise variable is taken to be & ): 
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Typical examples of the generalized forces in these equations 
are listed below: 


A) = the weighted life due to unit amplitude of bending 
th motion at the reference station, 


L 
7 AL LEG), (f<f at the reference Station), 


@.. = the weighted moment about the elastic axis due to 
Seat ame litude ef pitching motion at the reference 
stetion, 


- Lom) ally, 


and 
3) = the weighted moment about the elastic axis due to 
» unit amplitude of aileron motion at the reference 
station, 
ly ; 
Am é 
: ( eh |4y, 
where ( ),, denotes @ Woment about the elastic axis. The 


spanwise running lifts and moments, such as 44,@m,.).. and 
(4m,). etc. may be obtained by appropriste chordwise inte- 


grations of the pressures or of the velocity potentials which 
are calculated in advance by the box scheme. Finally, spanwise 
integrations may be carried out by numerical means. The quan- 


tities 7)... 7 are the elements of the flutter determinant 


\f 4-62 
22 
& wf 


Associated with this determinant are two unknowns (such as 
the velocity and the ae Sy At flutter these two unknowns 
( eigenvalues) are such that the value of the determinant is 


zero. 
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D.2 Plutter Equations for Plate-Type Structures 


For low-aspect-ratio surfaces or Plate-type structures, the 


assumption of beaa-rod type deformations is 


A new approach must be devised to treat these cases, 


ample of the modified procedure, consider a 


no longer justified. 
8 an ex- 
cantilever triangular 


wing which is divided into ten areas as shown in Fig. D.1. Each 
such area has a mass ™, concentrated at its Mente coast 


location ( Xf, ). If the structural influence coef 


x 
Fig 0.1 A Delta Wing Configuration with 
associated with these points are given, one may obtain by astrix 
iteration (Ref. 14) the first few free vibration frequencies and 


the relate orthogonal mode shapes. Once they are available, 
one may proceed to the flutter apnh ions of motion. If the Flut- 


ter mode is assumed to consist of the first 


icients 


Lumped Masses 


two free vibration 


modes A(t) /(%y¥) and BW) f # then the displacement of a 


point Xns fo) may be written as 


Bn lt) = MOL Onto) + BOY, bs fy) 


Eq. (D.2) 


For use in deriving Lagrange's equations, the total kinetic 


energy of the system is 
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4) 
oak 
KE « £ Zim é, 


Z ‘ty tee Bios ; 
= BAZ, Lhe ph] $B*Zirabf otal] + AB 2m, Moog] 
Eq. (D.3) 


Since the last term is zero by the Srehapent sty condition, the 
kinetic energy terms in Lagrange's equations are 


4 (sc bs At Zim, [fbco.tn)] 


a Q KE) = a + My ; Win ¥ 
4. 21H) ) = BD) lf bndes] 


Eqs. (D.4a-b) 
For simple harmonic motion, A (t) -wAlt), Blt) = -w* Be) . Because 


of the orthogonality of the modes, one may write for the poten- 
tial energy terms 


. 7 10 7 ‘ 
out = PAC) 24 1a h(t in) 


mre) » «3 Bie) Zim Lb Oth] - 
Fqe. (D.5a-b) 


Finally, the generalized forces @, J C, are: 
= []{lby+ (Pg) fora) aedy 
RY 


| LY ay 
s og if Mika. é ae jg (Ae f Hy, Nhlag)drdy #AQ 184, 
3 Eq. (D.6a) 
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se LN * yg Lh le pd dy 


, f ETL) 4) 4 4 4 BQ 5 
= [tops tt Lt hope se 1] Lhony » Ant 2a, 
wy Md Eq. (D.6b) 
where for inetance (AA) ay indicates the Pressure distribution 


due to the motion A) (xy) and the integrations are over the 


planform S , The task of evaluating these forces is discus: od 
in Appendix C. 


Equations (D.4a-b), (D.5a-b), and (D.6a-b) provide the 
elements in Lagrange's equations 


7 i) ae 9G 2 £48. a 


Equations gees form a set of two simulteneous homogeneous linear 
equations in the complex modal amplitude functions A and a 


ta 2 
(5 0), [flote)] hy Ray Na 
2 
2244 2 6 
Qe o DV rd bub), 
Eq. (D.8) 


The condition for flutter is thet the determinant of the above 
matrix should vanish. 
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